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ONCOMELANIA HUPENSIS (GASTROPODA: RISSOOIDEA), 
INTERMEDIATE HOST OF SCHISTOSOMA JAPONICUM IN CHINA: 
GENETICS AND MOLECULAR PHYLOGENY BASED ON 
AMPLIFIED FRAGMENT LENGTH POLYMORPHISMS 
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ABSTRACT 


Schistosomiasis japonica is one of the most serious parasitic diseases in China. Sub- 
species of the pomatiopsid snail species Oncomelania hupensis transmit the human blood 
fluke Schistosoma japonicum, the parasite causing the disease. In at least one study in- 
volving the subspecies Oncomelania hupensis robertsoni, there was no clear phyloge- 
netic concordance between the mitochondrial and nuclear data, and it was strongly 
recommended that future studies incorporate more data from nuclear loci in order to better 
understand phylogeography and host-parasite coevolution (Wilke et al., 2006). In this pa- 
per, we explore genetic diversity based AFLP data involving 25 populations throughout 
mainland China. AFLP is a DNA fingerprinting technique that detects the polymorphism of 
the whole genome without prior knowledge of the nucleotide sequence, with the high de- 
gree of reproducibility and small amount of template DNA needed. 

Our AFLP data show that: (1) the genetic diversity within an O. hupensis population 
increases gradually with snail sample size, and when the sample size is more than 30 
individuals, the genetic variation within one O. hupensis population trends to stabilization. 
(2) There is high intra-population genetic variation on the mainland of China, and these 
intra-population genetic variations from different areas differ significantly. (3) The consid- 
erable genetic differentiation occurs throughout China. (4) The genetic variation among 
populations of O. h. hupensis is higher than that of O. h. robertsoni. (5) The patterns of 
genetic differentiation are basically consistent with geographical distribution of snail popu- 
lations. Our results are similar to the results of allozymes and COI gene sequences, ex- 
cluding a snail population from Guangxi Province (Gx-1), indicating that there are three 
distinct subspecies in mainland China, namely O. h. hupensis, O. h. robertsoni, and O. h. 
tangi. However, our results do not support that the snail population from Guangxi Province 
belongs to O. h . hupensis. 

Key words: Schistosomiasis japonicum, Ocomelania hupensis, China, molecular genet- 
ics, phylogeny. 


INTRODUCTION 


Schistosomiasis japonica is one of the most 
serious parasitic diseases in China, with a 
documented history of over 2100 years. By 
the end of 2003, an estimated 843,011 people 
were infected with S. japonicum. Among them, 
there were 1,114 cases with acute schistoso- 
miasis and 24,441 cases suffering from ad- 
vanced schistosomiasis. The human 
population at risk of infection was estimated 


around 65 million, and there were 74,000 in- 
fected cattle. The intermediate snail host, On- 
comelania hupensis, infested an area of 3,787 
km? (Zhou et al., 20052). Basic epidemiologi- 
cal and malacological research were discour- 
aged in China as the theoretical goal of 
preventing schistosomiasis japonicum trans- 
mission seemed attainable by introduction of 
effective programs of chemotherapy (using 
Praziquantel) to eliminate adult worms in the 
human or livestock hosts, and proper sanita- 
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tion concurrent together with the indiscriminant 
application of molluscicides to exterminate the 
snails. Although such integrated programs 
have had dramatic success in some areas of 
China (Jiang et al., 2002), the sustained effort 
required for their maintenance in contiguous 
water bodies requires considerably more re- 
sources than are available in many affected 
areas (World Health Organization, 1993). Ad- 
ditionally, synthetic molluscicides can cause 
serious environmental damage to freshwater 
ecosystems. In China, water development 
schemes, such as the Three Gorges Dam on 
the Yangtze River and the south-to-north wa- 
ter transfer project, may encourage growth and 
spread of Oncomelania hupensis, creating a 
need for further basic malacological research. 

The subspecies and populations of On- 
comelania hupensis differ in shell size and 
structure, breeding environment, growth rates, 
population genetics, and potential for infection 
by S. japonicum (Yuan et al., 1984; Hope & 
McManus, 1994; Davis et al., 1995, 1999a, 
2006; Wilke et al., 2000; Zhou et al., 2005b). 
Based on congruence of allozymes and mito- 
chondrial cytochrome. oxidase | [COI] gene 
sequences (Davis et al., 1995, 1999a, b; Wilke 
et al., 2000; Shi et al., 2002), there are three 
distinct clades of Oncomelania hupensis in the 
mainland China, widely separated bio- 
geographically. These have been assigned to 
three subspecies. Oncomelania h. hupensis 
thrives throughout the Yangtze River drainage 
system below the Three Gorges to Shanghai, 
Jiangsu, and Zhejiang provinces (some popu- 
lations are dispersed at high altitudes). This 
subspecies has populations with smooth or 
ribbed shells (Davis et al., 2006). Oncomela- 
nia h. robertsoni is located in Sichuan and 
Yunnan provinces at high elevations (500 m, 
most 1,000 m above sea level) above the 
Three Gorges of the Yangtze River and has 
the ancestral (pleisiomorphic) small, smooth 
shell without a varix (thickening of the outer 
lip of the shell) (Davis, 1979; Davis et al., 
1995); O. h. tangi is found in Fujian Province 
along the southern coast of China, opposite 
Taiwan, separated from the Yangtze River by 
high mountain ranges. The shells are smooth 
but have a doubly thick varix, and their allom- 
etry differs from that of the other two subspe- 
cies in that their shells are proportionally much 
wider for their length. 

As was shown by previous allozyme stud- 
ies, the introduction of new molecular-based 
methodologies has the potential to greatly 
improve differentiating populations and higher 


taxa. In recent years, some universal PCR 
primers were used to amplify targets in nuclear 
and mitochondrial genomes across a broad 
range of taxa (Simon et al., 1994; Pedersen, 
1996), and the introduction of arbitrary primer 
mediated DNA fingerprinting (Williams et al., 
1990; Welsh & McClelland, 1990; Caetano- 
Anolles et al., 1991; Vos et al., 1995). The 
choice of these universal primer methods for 
species identification and elucidation of rela- 
tionships in tropical freshwater snails was dis- 
cussed by Jones et al. (1999). The Randomly 
Amplified Polymorphic DNAs (RAPDs) method 
was successfully utilized for species identifi- 
cation and clarification of genetic relationships 
among planorbid snails in Central and East 
Africa (Jones et al., 1999). In contrast, Ampli- 
fied Fragment Length polymorphisms (AFLPs) 
combine RAPDs and RFLPs to assay anony- 


mous sequences over the entire genome, but 


with the major advantage of being more ro- 
bust and reliable than other arbitrary primer- 
mediated DNA fingerprinting methods, 
because amplification is performed under 
more stringent reaction conditions (Vos et al., 
1995). The AFLP markers have been used to 
infer phylogenetic relationships of closely re- 
lated species based on measures of genetic 
distance (Heun et al., 1997; Semblat et al., 
1998). Based on analyses mtDNA, LSU rRNA, 
and AFLPs of O. hupensis robertsoni speci- 
mens, Wilke et al. (2006) found that there was 
no clear phylogenetic concordance between 
the mitochondrial and nuclear data in On- 
comelania hupensis robertsoni (in contrast 
with O. h. hupensis). They strongly recom- 
mended that future studies incorporate more 
data from nuclear loci in order to better under- 
stand phylogeography, population genetics, 
and host-parasite coevolution. 

The purposes of this paper are (1) to use 
AFLP methodology to assess the degree to 
which AFLPs derived relationships are con- 
gruent with the allozymes and COI results; (2) 
to determine the degree to which AFLPs data 
enables analysis of population structure. 


MATERIALS AND METHODS 
Locality Data 


Oncomelania is distributed in twelve prov- 
inces of mainland China. We collected 
samples of 25 snail populations from ten prov- 
inces. We did not find sufficient samples for 
studying in Shanghai and Guangdong prov- 
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FIG. 1. Localities studied. Populations pointed out by arrows were used for studies of intra-population 


variation. 


-` inces. A total of 760 uninfected individuals from 
25 O. hupensis populations were studied 
(Fig.1). Analyses were done on 40 individuals 
randomly sampled from the O. hupensis popu- 
lation from Hn-2; 30 individuals were analyzed 
from each of the other populations. Locality 
information is given in Table 1. 


Snail Selection 


Individuals (120) were collected at random 
from each O. hupensis population. The snails 
were isolated without food, in an antibiotic 
solution (ampicillin 100 ug/ml) and left over- 
night at 4°C to minimize prokaryotic growth 
prior to DNA extraction. Then the cercaria 
shedding test was used to examine whether 
they were infected with S. japonicum. Those 
positive snail individuals were eliminated. The 
remaining individuals were crushed and ex- 
amined microscopically for larval trematodes, 
then the negative individuals were selected for 
DNA isolation. 


DNA Extraction 


Genomic DNA was extracted from the foot 
of each snail individual using a modified CTAB 
protocol (Enjalbert et al., 2002). The DNA con- 
centrations were determined by O.D. measure- 
ment using a BioPhotometer (Eppendorf, 
Hamburg, Germany). All samples were then 


diluted to a concentration of 50ng/ul. After di- 
lution, additional measurements were per- 
formed to ensure the correct concentration. 


Population Used to Study Variation with In- 
creased Sample Size and Number of AFLP 
Molecular Loci 


The Hn-2 population was selected to ana- 
lyze how much variation increasing with sample 
size and the number of AFLP molecular loci. 


Populations Used to Study Intra-Population 
Genetic Variation 


Thirteen snail populations (Yn-1, Js-1, Sc-1, 
Fj-1, Gx-1, Js-2, Ah-1, Jx-1, Jx-2, Hn-1, Hn-2, 
Hb-1, Hb-2) were selected from the 25 snail 
populations to determine the amount of intra- 
population variation. Of these snail populations, 
ten populations (Hn-1, Hn-2, Js-1, Js-2, Gx-1, 
Hb-1, Hb-2, Ah-1, Jx-1 and Jx-2) from six prov- 
inces belong to O. h. hupensis; two popula- 
tions (Yn-1 and Sc-1) are O. h. robertsoni; one 
population (Fj-1) is O. h. tangi. 


DNA Pool Construction 


For each snail population, 30 ul from 30 snail 
individuals (1 ul per individual) from the same 
population were pooled together to assess the 
molecular genetics of O. hupensis. 


TABLE 1. Locality information for twenty five snail populations of Oncomelania hupensis studied. 


Locality code Shell sculpture Varix 


Hn-1 
Hn-2 
Hn-3 
Yn-1 
Js-1 
Js-2 
Js-3 
Sc-1 
Sc-2 
Sc-3 
Sc-4 
Sc-5 
Fj-1 
Gx-1 
Hb-1 
Hb-2 
Hb-3 
Hb-4 
Ah-1 
Ah-2 
Ah-3 
Ah-4 
Jx-1 
Jx-2 
Zj-1 


Ribbed 
Ribbed 
Ribbed 
Smooth 
Smooth 
Ribbed 
Ribbed 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Ribbed 
Ribbed 
Ribbed 
Ribbed 
Ribbed 
Smooth 
Ribbed 
Ribbed 
Ribbed 
Ribbed 
Ribbed 


Yes 
Yes 
Yes 

No 
Yes 
Yes 
Yes 

No 

No 

No 

No 

No 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 


Locality 


Hunan, Nanxian County, Babaigong Township, Xinshang Village 
Hunan, Yueyang City, Junshan District, Liulinzhou Township, Changjiang Village 
Hunan, Hua Rong County, Xing Fu Township, Dong Jia Village 
Yunnan, Dali City, Yingiao Township, Wuligiao Village 

Jiangsu, Yixing City, Hufu Township, Zhanggong Village 

Jiangsu, Dantu County, Shiye Township, Xintan Village 

Jiangsu, Jiang Ning District, Tong Jing Township, Xin Ji Zhou Village 
Sichuan, Pujiang County, Changqiu Township, Sanhe Village 
Sichuan, Dan Leng County, Yang Chang Township, Xu Ba Village 
Sichuan, Guang Han City, Shuang Quan Township, Long Quan Village 
Sichuan, Xi Chang City, Chuan Xing Township, Xin Long Village 
Sichuan, Pu Ge County, Te Zi Township, Chang Shou Village 

Fujian, Fuging City, Donghan Township, Donghan Village 

Guangxi, Yizhou City, Fulong Township, Jingkou Village 

Hubei, Wuhan City, Caidian District, Xiaosi Township, Wagou Village 
Hubei, Jiangling County, Zishi Township, Fanyuan Village 

Hubei, Han Chuan City, Liu Ge Township, Xi Si Tai Village 

Hubei, Yang Xin County, Mu Gang Township, Xin Hu Village 

Anhui, Tongling County, Laozhou Township, Guanghui Village 
Anhui, Tong Ling County, Shun An Township, Chang Shan Village 
Anhui, Cong Yang County, Feng Yi Township, Feng Yi Village 
Anhui, Gui Chi District, Mei Long Township, Guo Gang Village 
Jiangxi, Douchang County, Zhouxi Township, Chaipeng Village 
Jiangxi, Xingzi County, Xinchi Township, Zhuxi Village 

Zhejiang, Jiangshan County, Wangyao Township, Dahe Village 


Latitude 


29.1365°N 
29.3881"N 
29.3060°N 
25.6686°N 
31.2254°N 
32.2107°N 
31.7962 nN 
30.2167°N 
29:919 *N 
30.9952°N 
27.8670°N 
27.6129" IN 
25.4130°N 
24.5009°N 
30.4471°N 
30.1842°N 
30 756 53N 
29.7265 IN 
3059723 IN 
30:9538"N 
30.7106°N 
30.7673°N 
29.1822"IN 
29.3098"N 
28.6890°N 


Longitude 


112:3329°E 
112.8595" E 
112.7202" E 
100.1789: E 
11999871 E 
119.2635 È 
118.5447°E 
103:6370°E 
10335902 E 
104.4267°E 
102.329 E 
102:6659 E 
119.6140°E 
108.1668°E 
113.6669 E 
112.4084°E 
113.6144 E 
115.2476°E 
1T VETE 
117.9407°E 
117.4407°E 
117.6423°E 
116:39412*E 
116.0047^E 
11865974E 


OLE 


W X3POHZ 
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Amplified Fragment Length Polymorphism 
(AFLP) Reaction 


AFLP protocol was adapted from Mikkonen 
et al. (2005) with some modifications as fol- 
lows: two hundred nanograms of genomic DNA 
was digested with 8 units of each restriction 
Enzyme Pst/ and Msel (New England Biolabs), 
and then restriction fragments were ligated to 
the Pst! adapter and Msel adapter at 37°C for 
5 hin a restriction-ligation buffer [2.5ul 10 x Re- 
action buffer, 1 ul Pst// Msel adapter, 2 ul PstI 
/ Msel, 2.5 ul 10mM ATP, 3 unit of T4 DNA 
Ligase, 7 ul Ultrapure Water]. The final volume 
of the mixture was 20 yl. After a 1 : 20 dilution 
in ultrapure water, ligation products were used 
directly (no streptavidine purification) for the 
pre-amplification, using Pst and Mse primers 
(2 ng/uL) with no selective nucleotides. 

In the pre-amplification, 2 ul of ligated re- 
striction fragments was used as template. 1 ul 
unlabelled Pst/ -0 primer and Mse! -0 primer 
(MWG-Biotech AG) were added to a amplifi- 
cation buffer [1 ul dNTPs, 2.5 ul 10Reaction 
buffer, 1 unit of T4 DNA Ligase, 18 ul Ultrapure 
Water], giving a total volume of 25 ul. 
Preselective amplification was performed at 
94°C for 2 min, followed by 30 cycles of 94°C 
for 30 s, 56°C for 30 s and 72°C for 80 s and 
at last 72°C for 5 min. 

In the selective amplification, the products 
of the pre-amplification were diluted 1 : 20 in 
ultrapure water, and followed by a specific 
amplification with three selective nucleotides 
on both Mse! and Pst/ primers (Table 2). Se- 
lective DNA amplification was performed in a 
20 ul reaction mixture consisting of 2 ul of tem- 


plate DNA from the preselective PCR, 2 ul 10x 
Reaction buffer, 1 ul dNTP, 1 ul fluorescence 
labeled Pst/ primer (5 ng/ul), 1 ul fluorescence 
labeled Msel primer (30 ng/ul), 1 unit of T4 
DNA Ligase, and 12.5 ul Ultrapure Water. 
Selective amplification was done by incubat- 
ing at 94°C for 2 min, then one cycle of 94°C 
for 30 s, 65°C for 30 s and 72°C for 80 s. The 
annealing temperature was then reduced by 
0.7°C in each subsequent cycle until attaining 
56.6°C (12 cycles), followed by an extra 23 
cycles at a 94°C for 30 s, 55°C for 30, 72°C 
for 80 s, and at last 72°C for 5 min. 

PCR products were separated electrophoreti- 
cally using a 6% polyacrylamide gel and an 
automatic ABI 377 DNA Sequencer (GMI, Inc.). 
The resulting band patterns were analyzed vi- 
sually. AFLP fragments of strong intensity were 
scored as binary characters for each isolate. 


Data Analysis 


The AFLP markers, treated as alleles, were 
analyzed following the approach used for 
RAPD (Random Amplified Polymorphic DNA) 
markers (Luch, 1990; Vekemans et al., 2002). 
The percentage of polymorphic loci, Nei's gene 
diversity (Nei, 1973), and Shannon's informa- 
tion index (Lewontin et al., 1972) were com- 
puted using the software package PopGen32 
(version 1.32) (Yeh, 1999). Pairwise similari- 
ties (GS,,.) were calculated between snail 
populations using Dice's coefficient (Nei & Li, 
1979) by Simqual program of the software pack- 
ages NTSYpc (Rohlf, 1998), and the similarity 
coefficients (GS,...) were converted to genetic 
distance (Genetic distance D = 1-GS,, .).The 


TABLE 2. 12 primer combinations used for AFLP analysis, and their sequences. 


Primer Pair Pstl Primer Sequence 


Pre- PO  5-GACTGCGTACATGCA-3' 

AT P1  5-GACTGCGTACATGCAGAA-3' M7 
A8 P1  5-GACTGCGTACATGCAGAA-3' M8 
B4 P2 5-GACTGCGTACATGCAGAC-3' M4 
B6 P2 5-GACTGCGTACATGCAGAC-3' M6 
B7 P2 5-GACTGCGTACATGCAGAC-3' M7 
B8 P2 5-GACTGCGTACATGCAGAC-3' M8 
F5 P6 5-GACTGCGTACATGCAGTC-3' M5 
F7 P6 5-GACTGCGTACATGCAGTC-3' M7 
G5 P7  5-GACTGCGTACATGCAGTG-3' M5 
G7 P7  5-GACTGCGTACATGCAGTG-3' M7 
G8 P7 S-GACTGCGTACATGCAGTG-3’ M8 
H8 P8 5-GACTGCGTACATGCAGTT-3' M8 


Msel Primer Sequence 


MO 5-GATGAGTCCTGAGTAA-3' 
5-GATGAGTCCTGAGTAACTG-3' 
5-GATGAGTCCTGAGTAACTT-3 
5-GATGAGTCCTGAGTAACAT -3' 
9-GATGAGTCCTGAGTAACTC-3' 
5-GATGAGTCCTGAGTAACTOG-3' 
5-GATGAGTCCTGAGTAACTT-3' 
9-GATGAGTCCTGAGTAACTA-3' 
5-GATGAGTCCTGAGTAACTO-3' 
9-GATGAGTCCTGAGTAACTA-3' 
o9-GATGAGTCCTGAGTAACTG-3' 
o9-GATGAGTCCTGAGTAACTT-3' 
o9-GATGAGTCCTGAGTAACTT-3' 
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FIG. 2. Information on intra-population genetic diversity with snail sample size in- 
creasing. 
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FIG. 3. Information on intra-population genetic diversity with the number of AFLP loci 
increasing. 
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genetic diversities were compared among 
multiple snail populations with Kruskal-Wallis 
H nonparametric test, between subspecies of 
O. hupensis with Mann-Whitney U of the non- 
parametric test. GS,,, were compared be- 
tween subspecies of O. hupensis with 
bootstrap approach. The other data analysis 
and graphs drawing were performed using the 
software packages Microsoft Excel (Microsoft 
Corporation, USA) and SPSS (SPSS Inc., Chi- 
cago, USA). 


RESULTS 


Relationship between Stability of Information 
on Genetic Variation and Sample Size 


The AFLP method was used to analyze the 
40 individuals of O. hupensis from Hn-2 with 
12 primer pairs. 42—74 distinct fragments with 
a size in the range of 50-490 bp were yielded 
for each primer combination, and the percent- 
age of polymorphic loci was approximately 
82.9%-93.9%. Overall, the percentage of 
AFLP polymorphic loci, Nei's gene diversity, 
and Shannon's information index increased 
gradually with increasing sample size. There 
were three phases: (a) The rapidly changing 
phase: when sample size increased from 2 to 
7 snail individuals, the three indices above in- 
creased rapidly. (b) The slowly changing 
phase: when sample size increased from 7 to 
30 individuals, the values of the three indices 
increased slowly. (c) The stable phase: when 
sample size was more than 30 snails, the val- 


ues of these indices trended to stabilization 
(Fig. 2). The snail individuals in the order on 
the x-axis were randomized. 


Relationship Between Stability of Information 
on Genetic Variation and Number of AFLP 
Molecular Loci 


Thirty out of 40 snails from Hn-2 above were 
selected at random to analyze how much 
variation increased with the number of AFLP 
loci. With the increase in the number of AFLP 
loci, the change of the percentage of polymor- 
phic loci and the indices of genetic diversity 
was divided approximately into three phases 
as above: (a) The phase of rapid change: when 
the number of AFLP loci increased from 2 to 
128, the values of the percentage of polymor- 
phic loci and the indices of genetic diversity 
fluctuated rapidly. (b) The phase of slow 
change: when it increased from 128 to 338, 
the values of these indices fluctuated slowly. 
(c) The stable phase: when the number of 
AFLP loci was more than 338, the values of 
these indices trended to stabilization (Fig. 3). 
The loci in the order on the x-axis were ran- 
domized. 


Intra-Population Genetic Diversity 


Thirteen populations (Fig. 1) were studied 
by AFLP method with eight primer pairs 
(namely A7, A8, B4, B6, B7, G5, G7, and G8). 
Within single locality population, the number 
of fragments amplified by AFLP method 
ranged from 403 to 472 with the number of 


TABLE 3. Genetic diversities within 13 populations of Oncomelania hupensis. 


Population No. No. of Percentage of Nei's gene Shannon's 
code loci polymorphic loci polymorphic loci (96) diversity information index 
Yn-1 436 379 86.87 0.322 + 0.0026 0.476 + 0.0035 
Js-1 403 314 77.92 0.270 + 0.0030 0.404 + 0.0041 
Sc-1 424 314 74.03 0:265 0:0031 0.421 + 0.0043 
Fj-1 430 304 70:61 0.253 + 0.0029 0.377 + 0.0042 
Gx-1 405 226 55.80 0.901 + 0.0031 0.287 + 0.0045 
Js-2 467 418 89.51 0.333 0.0028 0.493 + 0.0031 
Ah-1 426 370 86.85 0.322 + 0.0025 0.478 + 0.0035 
Jx-1 434 378 87.10 0.313 + 0.0026 0.465 + 0.0035 
Jx-2 472 440 98.22 0.345 + 0.0022 0.510 + 0.0029 
Hn-1 459 403 87:77 0.339 + 0.0025 0.500 + 0.0034 
Hn-2 453 391 86.31 0.276 + 0.0027 0.420 + 0.0037 
Hb-1 433 9T 85.88 0.330 + 0.0026 0.485 + 0.0035 
Hb-2 415 343 82.55 0.302 + 0.0029 0.450 + 0.0040 
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polymorphic loci from 226 to 440; a percent- 
age of polymorphic loci ranged from 55.80% 
to 93.22%, and the Nei’s gene diversity from 
0.191 to 0.345; the Shannon’s information in- 
dex ranged from 0.287 to 0.510. The genetic 
diversities within Oncomelania populations 
from different areas differed significantly (Nei’s 
gene diversity y? = 2500.56, P < 0.01; 
Shannon’s information index y? = 2501.20, P 
< 0.01). Of thirteen snail populations, the ge- 
netic diversity within Gx-1 was the lowest, and 
the percentage of polymorphic loci, Nei’s gene 
diversity, Shannon’s information index and 
average coefficient of similarity were 55.80%, 
0.191, 0.287 and 0.904 respectively; whereas 
the percentage of polymorphic loci, Nei’s gene 
diversity and Shannon's information index of 
Jx-2 were the most and were 93.22%, 0.345 
and 0.510 respectively (Table 3). 


Inter-Population Genetic Variation 


DNA pools of 25 five snail populations were 
amplified using the AFLP method with eight 
primer pairs (A7, A8, B4, B6, B7, G5, G7, and 
G8). These distinct amplified fragments with 
a size in the range of 50 to 325 bp were se- 
lected to analyze the genetic variation among 
Snail populations. 

Among 25 snail populations, Nei’s gene di- 
versity and Shannon’s information index were 
0.280 and 0.426 respectively; genetic distance 
D from 0.169 to 0.306, and Nei’s unbiased 
genetic distance from 0.225 to 0.452. Among 
six snail populations of O. h. robertsoni, ge- 
netic distance D was from 0.170 to 0.235, and 
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Nei's unbiased genetic distance was from 
0.225 to 0.334. Among 18 snail populations of 
O. h. hupensis, genetic distance D was in the 
range of 0.176 to 0.289, and Nei's unbiased 
genetic distance was in the range of 0.229 to 
0.413. Whether Nei's gene diversity, 
Shannon’s information index, genetic distance 
D, or Nei's unbiased genetic distance, among- 
snail population distance of O. h. hupensis was 
greater than that between O. h. robertsnoi 
populations (P « 0.01). Of 18 snail populations 
of O. h. hupensis, 15 snail populations had 
ribbed-shells, and three snail populations had 
smooth shells. Nei's gene diversity and 
Shannon's information index among ribbed- 
shell populations of O. h. hupensis were higher 
than that among smooth-shelled populations 
of O. h. hupensis (P « 0.01), but genetic dis- 
tance D and Nei's unbiased genetic distance 
among ribbed-shelled populations of O. h. 
hupensis were lower than that among smooth- 
shelled populations of O. h. hupensis (P « 
0.01). If the smooth-shelled populations of O. 
h. hupensis were not included a snail popula- 
tion from Guangxi province, Genetic distance 
D and Nei's unbiased genetic distance among 
smooth-shelled populations of O. h. hupensis 
were not different from that among ribbed- 
shelled populations of O. h. hupensis (P > 
0.05). Nei's gene diversity and Shannon's in- 
formation index among populations of O. h. 
robertsoni were not different from that among 
smooth-shell populations of O. h. hupensis (P 
> 0.05), but genetic distance D and Nei’s un- 
biased genetic distance among smooth-shell 
populations of O. h. hupensis were higher than 


TABLE 4. Genetic variation among the populations of Oncomelania hupensis. 


O. h. hupensis Total 


Smooth Total 


0.209 € 0.0035 0.262+0.0027 0.280 + 0.0026 


0.300 x 0.0051 0.399+0.0037 0.426 + 0.0035 


O. h. robertsoni 
Ribbed 
Nei's genetic 0.201 + 0.0030 0.253 + 0.0028 
diversity 
Shannon's — 0.300 + 0.0034 0.384 + 0.0039 
information index 
Nei's unbiased 0.278 + 0.0088 0.310 x 0.0035 
genetic distance 
Neil's unbiased 0.758 + 0.0067 0.729 + 0.0026 
genetic identity 
Sbe 0.797 + 0.0055 0.765 + 0.0022 
Genetic distance D 0.203 + 0.0055 0.235 + 0.0022 


0.380 + 0.0150 


0.684 + 0.0101 


0.722 + 0.0064 
0.278 + 0.0064 


0.328 + 0.0032 


0.722 + 0.0023 


0.759 + 0.0020 
0.241 + 0.0020 


0.328 + 0.0039 


0.722 + 0.0028 


0.761 + 0.0024 
0.239 + 0.0024 
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FIG. 4. Complete-link trees depicting phylogenetic relationships among twenty five snail populations 
of O. hupensis using (a) Genetic distance D matrix and (b) Nei's unbiased genetic distance matrix. 


that among populations of O. h. robertsoni (P 
« 0.01). Whether Nei's gene diversity, 
shannon's information index, Genetic distance 
D or Nei's unbiased genetic distance among 
ribbed-shell populations of O. h. hupensis was 
larger than that among snail populations of O. 
h. robertsoni (P « 0.01) (Table 4). 


Snail AFLP Phylogeny 


Analyses of Genetic distance D and Nei's 
unbiased genetic distance matrixes using the 
Complete-link method, Neighbor-joining (NJ) 
method, Single-link method and Unweighted 
pair-group method with arithmetic means 


3/6 


0.00 0.04 
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Fj-1 


Hb-4 


Ah-4 


0.08 


Genetic distance D 


b) 
Fj-1 
Gx-1i 
Yn-1 
Sc-2 
Sc-1 
Sc-3 
Sc-4 
Sc-5 
Hn-3 
Hn-2 
Hn-1 
Jx-1 
Jx-2 
4s-2 
Hb-1 
Hb-3 
Js-3 
Js-1 
Hb-4 
Zj-1 
Hb-2 
Ah-3 
Ah-4 
Ah-1 
Ah-2 
0.00 9.06 0.12 0.18 0.23 


Nei's unbiased genetic distance 


FIG. 5. Neighbor-joining trees depicting phylogen 


etic relationships among twenty five snail populations 


of O. hupensis using (a) Genetic distance D matrix and (b) Nei's unbiased genetic distance matrix. 


(UPGMA) method respectively were shown in 
Figs. 4-7. For the complete-link method and 
NJ method, whether Genetic distance D ma- 
trix or Nei's unbiased genetic distance matrix 
was analyzed, the two Complete-link trees and 
the two NJ trees were similar to each other, 
and especially the NJ phylogram using Ge- 


netic distance D matrix was very similar to the 
NJ tree using Nei's unbiased genetic distance 
matrix. These four trees were all rooted with 
Fj-1 and Gx-1 populations and showed two 
big branches. One branch included six 
smooth-shelled populations without varix (Yn- 
1, Sc-1-Sc-5) from Yunnan and Sichuan Prov- 
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FIG. 6. UPGMA trees depicting phylogenetic relationships among twenty five snail populations of O. 
hupensis using (a) Genetic distance D matrix and (b) Nei’s unbiased genetic distance matrix. 


inces, the other branch consisted of other 17 
snail populations from throughout the Yangtze 
River drainage below the Three Gorges to 
Jiangsu and Zhejiang provinces. The two 
UPGMA trees also were both similar to the 
four trees above, and were both rooted with 
Fj-1 and Gx-1 populations and showed two 


big branches. Unexpectedly, three ribbed- 
shelled populations (Hn-1, Hn-2 and Hn-3) 
from Hunan Province fell within the O. 
hupensis robertsoni subspecies (which is 
found in Yunnan and Sichuan provinces) clade 
in the UPGMA tree using Genetic distance D 
matrix, but the genetic distance D between 
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Fj-1 
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Ah-3 
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FIG. 7. Single-link trees depicting phylogenetic relationships among twenty five snail populations of 
O. hupensis using (a) Genetic distance D matrix and (b) Nei’s unbiased genetic distance matrix. 


the snail populations from Hunan Province and 
the snail populations from Yunnan and Sichuan 
provinces was relatively big (0.220 to 0.298). 
The two single-link trees were different from 
six other trees above, and were both rooted 
with only Fj-1 populations, with Gx-1 popula- 


tion divided into an independent branch. The 
other 23 snail populations also were divided 
similarly into two big branches except the Js- 
1 population in the tree using genetic distance 
D matrix, which was divided into an indepen- 
dent branch. 
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DISCUSSION 


Appropriate Sample Size and Number of AFLP 
Markers 


Due to intra-population genetic variation, one 
or several snail individuals are insufficient to 
analyze population variation. The snail popu- 
lation (Hn-2) selected to examine how much 
variation increased with sample size, and the 
number of AFLP molecular loci showed that 
the percentage of polymorphic loci, Nei's gene 
diversity, and Shannon's information index all 
increased gradually with increasing sample 
size up to a sample size of 30 or more indi- 
viduals, When the sample size was less than 
seven individuals, the values of these indices 
are very unstable. Our result is consistent with 
the sample size calculated based on the for- 
mula of Brown & Moran (1981) in a random 
mating population. So we sampled randomly 
30 snail individuals not infected with S. 
japonicum from a snail population to represent 
this population in this paper. We recognize the 
fact that if additional populations were added, 
variation would continue to increase. The prac- 
tical solution, universally recognized, is to use 
between 15 to 25 individuals for population 
genetic studies. 

Our AFLP analysis data still showed, on ba- 
sis of 30 snails sample from Hn-2 population, 
that the percentage of polymorphic loci, Nei's 
gene diversity, and Shannon's information in- 
dex change rather acutely when the number 
of loci is less than 128, while the values of 
theses indices trend to stabilization when the 
number of loci is more than 338. These indi- 
cate that the appropriate number of AFLP 
markers for studying the genetic variations of 
O. h. by AFLP method is optimally not less 
than 338 for each Oncomelania population, 
whereas the result is very unreliable when the 
number of AFLP markers is less than 128. In 
our study, the total eight primer pairs were used 
to analyze genetic diversity of O. hupensis, 
and the number of AFLP markers amplified 
was in range of 403—467 bands for the 13 
Oncomelania populations. So it was enough 
and reasonable to use only eight primer pairs 
for studying and comparing the genetic varia- 
tion of difference Oncomelania populations. 


Intra-Population Genetic Diversity 


In a previous allozyme study of O. hupensis, 
Davis et al. (1995) found that 14 snail popula- 
tions from Mainland China had few alleles per 
locus (1.0 to 1.6), a low mean heterozygosity 


(0.008 to 0.036), and low percentage polymor- 
phic loci (4 to 28). In another previous COI gene 
sequencing study of O. hupensis, Wilke et al. 
(2000) also found relatively little genetic vari- 
ability within the in-group individuals, and the 
nucleotide diversity ranged from O to 0.0196. 
Compared with these relatively conservative 
allozyme and COI marker, our AFLP data 
showed that there was higher intra-population 
genetic variation in mainland China — the per- 
centage of polymorphic loci ranged from 
55.80% to 93.22%, the Nei's gene diversity 
from 0.191 to 0.345, and the Shannon’s infor- 
mation index from 0.287 to 0.510 — further- 
more, these intra-population genetic variations 
from different areas differed significantly. These 
substantial differences were present among not 
only the three subspecies, but also different 
populations of the same subspecies, especially 
O. h. hupensis. The results found in our study 
were similar to those by Wilke et al. (2006), 
and they found each snail individual had a 
unique AFLP fingerprint in the preliminary AFLP 
analysis of 24 specimens. So our results might 
indicate that AFLP approach can examine more 
genetic diversity of O. hupensis than allozyme 
and COI gene sequence. 

Oncomelania hupensis tangi has been eradi- 
cated except for two known populations, and 
the parasite presumably is extinct. Our AFLP 
data showed interestingly that the genetic 
variation within Fj-1 from O. h. tangi as well as 
Gx-1 population was very low. Schistosomia- 
Sis Japonicum in Guangxi Province has been 
exterminated, except for intermediate snail 
hosts restricted in narrow areas. This might be 
mostly due to the fact that these snail popula- 
tions are localized relatively narrow areas, and 
the number of individuals in these snail popu- 
lations is very limited, so the genetic drift can 
easily happen in these snail populations, and 
the results of the genetic drift will lead to far- 
ther decrease of their genetic diversity. 


Inter-Population Genetic Diversity and Phylog- 
eny 


Our AFLP analysis data from 25 snail popu- 
lations show that considerable genetic differ- 
entiation has been occurred throughout China. 
Genetic distance D was from 0.169 to 0.306, 
and Nei's unbiased genetic distance from 
0.225 to 0.452. These two indices of genetic 
distance with the four clustering methods were 
used to reconstruct the phylograms of O. 
hupensis from mainland China, and most re- 
sults of the eight phylogenetic trees were simi- 
lar, and patterns of genetic differentiation were 
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in general correlated with geographic distribu- 
tion. These results were similar to the results 
of allozymes and COI gene sequences (Davis 
et al., 1995, 1999a, b; Wilke et al., 2000; Shi 
et al., 2002) excluding a snail population from 
Guangxi Province (Gx-1). These patterns thus 
suggested that there were three distinct clades 
of O. hupensis on mainland China. 

The six snail populations from Yunnan and 
Sichuan provinces, namely O. h. robertsoni, 
were less, and they always fell within a clade 
in eight phylograms. This supported that snails 
from Yunnan and Sichuan provinces belong to 
the same subspecies, O. hupensis robertsoni. 

Excluding in two single-link trees, Fj-1 and 
Gx-1, populations always fell within a clade in 
other six trees. This was dissimilar to the re- 
sult of allozyme studies (Davis et al., 1995), 
namely our AFLP analysis data did not sup- 
port that the Gx-1 population belongs to O. 
hupensis hupensis. In fact, the Gx-1 popula- 
tion is far from the other populations of O. h. 
hupensis, and no snail populations have also 
been found in the areas between them. The 
genetic distance D between Gx-1 and Fj-1 snail 
populations was 0.240, and Nei's unbiased 
genetic distance was 0.342. There are some 
morphological similarities between Fj-1 and 
Gx-1 populations: the shells of these snails are 
smooth with a varix, and proportionally wider 
for their length. But the geographical distance 
between the two snail populations is far, and 
the shells of Gx-1 snails are thin with weak varix 
compared with that of Fj-1 population. The 
genetic distance between Fj-1 and Gx-1 popu- 
lations is an average level of that among 25 
populations. So, should the snails from 
Guangxi Province be an independent subspe- 
cies (i.e., O. hupensis guangxiensis) accord- 
ing to the classification of Liu et al. (1981), or 
be incorporated with snails from Fujian Prov- 
ince into a subspecies? This need further study. 

Among the other 17 snail populations, the 
Js-1 population was an independent clade in 
the single-link tree using genetic distance D 
matrix, and the three ribbed populations (Hn- 
1, Hn-2 and Hn-3) fell within the O. hupensis 
robertsoni subspecies clade in the UPGMA 
tree using genetic distance D. matrix. How- 
ever, these snail populations always fell within 
a clade in the other six trees. These results 
were similar to the results of allozyme and COI 
gene sequence, and showed the 17 popula- 
tions are a subspecies, O. hupensis hupensis. 

Our results still show the genetic variation 
between snail populations of O. h. hupensis is 
greater than that between populations of O. h. 


robertsoni. This might indicate there is bigger 
divergence among O. h. hupensis, compared 
with O. h. robertsoni. However, Wilke et al. 
(2006) found that in their study of 66 individu- 
als from 13 populations of O. h. robertsoni, that 
genetic divergence, using combined COI/LSU 
rRNA, was an order of magnitude greater than 
among populations of O. h. hupensis (0.042 
and 0.085 for different lineages of O. h. 
robertsoni; 0.0097 for other O. h. subspecies. 
But when AFLP was used on 24 individuals, 
each had a unique fingerprint. 


ACKNOWLEDGEMENTS 


This work was supported by the National 
Natural Science Foundation of China (No. 
30671799), in part by the S&T programs dur- 
ing the Tenth Five-year Plan (No. 
2004BA718B04) and in part by the National 
Natural Science Foundation of China 
(No.3059374). We are very grateful to Prof. 
George M. Davis for reviewing and modifying 
a draft of this paper, offering sound advice, 
and providing some literatures. We, the au- 
thors take responsibility for the final result. 


LITERATURE CITED 


BROWN, A. H. D. & G. F. MORAN, 1981, 
Isozymes and the genetic resources of forest 
trees. Pp. 1-10, in: M. T. CONKLE, Isozymes of 
North American forest trees and forest insects. 
Pacific South West Forest and Range Experi- 
ment Station General Technical Report, 48. 
U.S. Forest Service, Department of Agriculture. 


64 pp. 

CAETANO-ANOLLES, G., G. J. BASSAM & P. 
M. GRESSHOF, 1991, High resolution DNA 
amplification fingerprinting using very short 
arbitrary oligonucleotide primers. Biotechnol- 
ogy, 9: 553-556. 

DAVIS, G. M., 1979, The origin and evolution of 
the gastropod family Pomatiopsidae, with em- 
phasis on the Mekong River Triculinae. Mono- 
graph of the Academy of Natural Sciences of 
Philadelphia, 20:1—120. 

DAVIS, G. M., T. WILKE, W. P. WU & X. J. XU, 
2006, Ecogenetics of shell sculpture in On- 
comelania (Gastropoda) in canals of Hubei, 
China, and relevance for Schistosome trans- 
mission. Malacologia, 48(1—2): 253-264. 

DAVIS, G. M., T. WILKE, Y. ZHANG, X. XU, C. P. 
QIU, C-NLSPOLSIKY, D. C. Olu, Y. S. Ll, M. Y. 
XIA & Z. FENG, 1999a, Snail-Schistosoma, Para- 
gonimus interactions in China: population ecol- 
ogy, genetic diversity, coevolution and emerging 
diseases. Malacologia, 41(2): 355-377. 

DAVIS, G. M., Y. ZHANG, Y. H. GUO & C. M. 
SPOLSKY, 1995, Population genetics and sys- 


AFLP GENETICS OF ONCOMELANIA HUPENSIS 381 


tematic status of Oncomelania hupensis (Gas- 
tropoda: Pomatiopsidae) throughout China. 
Malacologia, 37(1): 133-156. 

DAVIS, G. M., Y. ZHANG, X. J. XU & X. X. Yang, 
1999b, Allozyme analyses test the taxonomic 
relevance of ribbing in Chinese Oncomelania 
(Gastropoda: Rissoacea: Pomatiopsidae). 
Malacologia, 41(1): 297-317. 

ENJALBERT, J., X. DUAN, T. GIRAUD, D. 
VAUTRIN, C. DE VALLAVIEILLE-POPE & M. 
SOLIGNAC, 2002, Isolation of twelve micro- 
satellite loci, using an enrichment protocol, in 
the phytopathogenic fungus Puccinia striiformis 
f. sp. tritici. Molecular Ecology Notes, 2: 563- 
505. 

HEUN, M., R. SCHAFER-PREGL, D. KLAWAN, 
R. CASTAGNA, M. ACCERBI, B. BORGHI & 
F. SALAMINI, 1997, Site of einkorn wheat do- 
mestication identified by DNA fingerprinting. 
Science, 278: 1312-1314. 

HOPE, M. & D. P. MCMANUS, 1994, Genetic 
variation in geographically isolated populations 
and subspecies of Oncomelania hupensis de- 
termined by a PCR-based RFLP method. Acta 
Tropica, 57: 75-82. 

JIANG, Q. W., L. Y. WANG, J. G GUO, M. G. 
CHEN, X. N. ZHOU & D. ENGELS, 2002, Mor- 
bidity control of schistosomiasis in China. Acta 
Tropica, 82: 115-125. 

JONES, C. &, L. R. NOBLE, J. OUMA, Fk C. 
KARIUKI, R. MIMPFOUDI, D. S. BROWN & D. 
ROLLINSON, 1999, Molecular identification of 
schistosome intermediate hosts: case studies 
of Bulinus forskalii group species (Gastropoda: 
Planorbidae) from Central and East Africa. Bio- 
logical Journal of the Linnean Society, 68: 215- 
240. 

LEWONTIN, R. C., 1972, The apportionment of 
human diversity. Evolutionary Biology, 6: 381- 
398. 

LIU, Y. Y., T. K. LOU, Y. X. WANG & W. Z. ZHANG, 
1981, Subspecific differentiation ‘of 
Oncomelaniid snails. Acta Zootaxonmica 
Sinica, 6: 253-266. 

LUCH, M., 1990, The similarity index and DNA 
finger printing. Molecular Biology and Evolu- 
tion, 7: 478-484. 

MIKKONEN, T., J. M. K. KOORT, K. J. 
BJORKROTH & A. SUKURA, 2005, Testing of 
amplified fragment length polymorphism 
(AFLP) technique as a tool for molecular epi- 
demiology. Veterinary Parasitology, 132: 19-22. 

NEI, M., 1973. Analysis of gene diversity in sub- 
divided populations. Proceedings of the Na- 
tional Academy of Sciences of the United 
States of America, 70: 3321-3323. 

NEI, M. & W. H. LI, 1979, Mathematical model 
for studying genetic variation in terms of re- 
striction endonucleases. Proceedings of the 
National Academy of Sciences of the United 
States of America, 76: 5269-5273. 

PEDERSEN, B. V., 1996, A phylogenetic analy- 
sis of cuckoo bumblebees (Psithyrus, 
Lepeletier) and bumblebees (Bombus, 
Latreille) inferred from sequences of the mito- 
chondrial gene cytochrome oxidase |. Molecu- 
lar Phylogenetics and Evolution, 5: 289-297. 


ROHLF, F. J., 1998, NTSYS-pc, numerical tax- 
onomy and multivariate analysis system, ver 
2.02. Exeter Software, Setauket, N.Y. 

SEMBLAT, J. P., E. WAJNBERG, A. DALMASSO, 
P. ABAD & P. CASTAGNONE-SERENO, 1998, 
High-resolution DNA fingerprinting of parthe- 
nogenetic root-knot nematodes using AFLP 
analysis. Molecular Ecology, 7: 119-125. 

SHI, C. H., T. WILKE, G. M. DAVIS, M. Y. XIA & 
C. P. QIU, 2002, Population genetics, micro- 
phylogeography, ecology and susceptibility to 
schistosome infection of Chinese Oncomela- 
nia hupensis hupensis (Gastropoda: 
Rissooidea: Pomatiopsidae) in the Miao River 
System. Malacologia, 44(2): 333-347. 

SIMON, C., F. FRATI, A. BECKENBACH, B. 
CRESPI, H. LIU & P. FLOOK, 1994, Evolution, 
weighting, and phylogenetic utility of mitochon- 
drial gene sequences and a compilation of 
conserved polymerase chain reaction primers. 
Annuals of the Entomological Society of 
America, 87: 651—701. 

VEKEMANS, X., T. BEAUWENS, M. LEMAIRE 
& |. ROLDAN-RUIZ, 2002, Data from amplified 
fragment length polymorphism (AFLP) mark- 
ers show indication of size homoplasy and of a 
relationship between degree of homoplasy and 
fragment size. Molecular Ecology, 11: 139-151. 

MOS, P, R. HOGERS, M. BLEBRER, M. 
REIJANS, T. VAN DE LEE, M. HORNES, A. 
PRIJTERS, J"POT, J! PELEMAM, M. KUIPER 
& M. ZABEAU, 1995, AFLP: a new technique 
for DNA fingerprinting. Nucleic Acids Research, 
11: 4407-4414. 

WELSH, J. & M. MCCLELLAND, 1990, Finger- 
printing genomes using PCR with arbitrary 
primers. Nucleic Acids Research, 18: 7213- 
"165. 

WILKE, T.. G M. DAVIS, C. E. CHEM X. WN. 
ZHOU, X. P. ZENG, Y. ZHANG & €. M. 
SPOLSKY, 2000, Oncomelania hupensis (Gas- 
tropoda: Rissooidea) in eastern China: molecu- 
lar phylogeny, population structure, and 
ecology. Acta Tropica, 77: 215-227. 

WILKE, T., G. M. DAVIS, D. C. Oli &Is..C. 
SPEAK, 2006, Extreme mitochondrial se- 
quence diversity in the intermediate schisto- 
somiasis host Oncomelania hupensis 
robertsoni: another case of ancestral polymor- 
phism? Malacologia, 48(1—2): 143-157. 

WILLIAMS, J. G. K., A. R. KUBELIK, K. J. LIVAK, 
J. A. RAFOLSKI & S. V. TINGEY, 1990, DNA 
polymorphisms amplified by arbitrary primers 
are useful as genetic markers. Nucleic Acids 
Research, 18: 6531-6535. 

WORLD HEALTH ORGANISATION, 1993, The 
control of schistosomiasis. WHO Technical 
Report Series 830. Geneva, Switzerland. 

YER, F. C., R. YANG & T. BOYLE, 1989, 
PopGene version 1.32. Microsoft Windows- 
based Freeware for population genetic analy- 
sis. University of Alberta, Edmonton, Alberta. 

YUAN, H. C., E. S. UPATHAM, M. KRUATRA- 
CHUE & V. KHUNBORIVAN, 1984, Suscepti- 
bility of snail vectors to Oriental anthropophilic 
Schistosoma. Southeast Asian Journal Tropi- 
cal Medicine and Public Health, 15: 86-94. 


382 ZHOU ET AL. 


ZHOU, X. N., L. Y. WANG, M. G. CHEN, X. H. 
WU, W. Q. JIANG, Y. X. CHEN, J. ZHENG & J. 
UTZINGER, 2005a, The public health signifi- 
cance and control of schistosomiasis in China 
— then and now. Acta Tropica, 96: 97-105. 

ZHOU, Y. B., Q. W. JIANG, G. M. ZHAO & J. G. 
WEI, 2005b, Analysis of morphological varia- 
tion within Oncomelania hupensis population. 
Chinese Journal of Zoology, 40: 90-97. 


Revised ms. accepted 25 February 2007 


